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ABSTRACT 

The possibility of velocity shear-induced linear transformations of different 
magnetohydrodynamic waves in the solar wind is studied both analytically and 
numerically. A quantitative analysis of the wave transformation processes for all 
possible plasma-/? regimes is performed. By applying the obtained criteria for 
effective wave coupling to the solar wind parameters, we show that velocity shear- 
induced linear transformations of Alfven waves into magneto-acoustic waves could 
effectively take place for the relatively low-frequency Alfven waves in the energy 
containing interval. The obtained results are in a good qualitative agreement 
with the observed features of density perturbations in the solar wind. 



Subject headings: solar wind: general, wave coupling 
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Introduction 



In spite of several decades of intense studies, the problem of the connection between 
the physics of the inner parts of the solar atmosphere and the solar wind remains incom- 
pletely understood. There is a general consensus that the mechanisms of coronal heating 
are somehow related to the mechanisms responsible for the acceleration of the solar wind. 
The so-called "basal" coronal heating (at r < 1.5 Rq) is usually attributed to a mixture 
of processes, e.g. wave dissipation (by phase mixi ng and /or resonant absorption) , magnetic 
reconnection, turbulence and plasma instabilities ( Cranmerll2002l ; ICranmer 112004 ) . In 'open' 
magnetic flux tubes which are feeding the fast solar wind at r > 2 Rq where the solar 
plasma is largely collision-less, however, additional heating seems to take place. This con- 
clusion follows from a number of bona fide observational signatures, viz. (a) low electron 
temperatures (T e < 1.5 x 10 6 K) in coronal holes, (b) in situ observations of the tempera- 
ture anisotropy T p > T e at 1 AU, (c) low radial gradients of these temperatures throughout 
the solar wind. These signatures unequivocally indicate at a gradua l, temporally and spa- 
tially extended addition of energy (ICranmer &: van Ballegooijenll2003l ) to the solar wind (the 
so-called "extended coronal heating" ) at a wide range of heliocentric distances. 

Presumably, propagating perturbations - magnetohydrodynamic (MHD) waves, vor- 
tices, shocks, turbulent eddies - are the best candidates for the transmission of energy from 
the inner to the outer parts of the solar atmosphere and further into the solar wind. In 
particular, Alfven waves are believed to play a significant role in the coronal h eating and 
the subsequent acceleration of the solar wind (ICranmer &: van Ballegooijenll2005l ). They are 
also considered as important diagnostic tools of the various physical processes occurring in 
solar plasmas. 

The presence of these efficient "energy-transmitters" is currently quite convincingly 
established. Even though previous observational studies only emphasized the appearance 
of Alfven waves in the solar atmosphere, in the modern era of satellite-based solar studies 
it became increasingly clear that the solar atmosphere hosts all three basic MHD wave 
modes, i.e. Alfven waves and both slow and fast magneto-acoustic waves. Many different 
sets of observations contributed to this ever growing evidence: in coronal loops the evidence 
came from the Extreme Ultra Violet (EUV) data from the TRACE satellite; in coronal 
plumes (slow magneto-acoustic waves) from the EIT instrument on-board SOHO; in the solar 
wind (propagat ing Alfven waves) from in-situ Helios and Ulysses spacecraft measurements 



rtCranmer 



20021 ). Therefore, both the observational data and the current theoretical studies 



are strongly in favor of the presence of a causal bond between the heating of the solar corona 
and the acceleration of the solar wind and some, probably significant, role of the MHD waves 
in the functioning of this bond. However, it is still unknown up to what extent the fluxes 
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of the solar wind's mass, momentum, and energy are driven by "classic", Parker-type gas 
pressure gradients, and what is the relative share of wave-flow (wave-particle) interactions 
and/or turbulence (jCranmer 112004 ). 



It is widely believed that one of the main mechanisms responsible for the heating of 
the solar wind, is the turbulent cascade of Alfvenic fluctuations. The advantage of this 
heating model is t hat it can also explain the above-mentio ned relatively high proton tem- 
peratures at 1 AU ([Goldstein. Roberts fc Matthaeus Ill995l ). Remarkably, it was also shown 



([Goldstein fc Roberts 



19991 ) that the turbulent cascade in the solar wind seems to evolve 



most rapidly in areas with a formidable velocity shear. At the s ame time, the measurements 

of sig nificant density fluctuations (Sp/p ~ 0.1) in the solar wind (IGoldstein. Roberts fc Matthaeus 
19951 ) are usually interpreted as evidence for the presence of magneto-acoustic waves not only 
in the corona but also in the solar wind. On the other hand, both slow and fast magneto- 
acoustic modes ar e expected to be strongly damped by Landau damping for solar wind 
parameter values ([Barnes Ill979l ). Therefore, there is a very scarce, if any, chance for their 
undamped propagation from the solar corona up to the outer wind regions. This logically 
leads to the following presumable solution of this puzzle: the magneto-acoustic waves in the 
solar wind should have a local origin. 

It has been suggested that three- an d /or four- wave resonant processes may be of a con- 
sider a ble importance i n the solar wind (jLacombe fc Mangeney I Il980l ; iBhattacharjee fc Ng 
200ll ; IChandran 1 120051 ). However, after studying the nonlinear interaction of oblique fast 
magneto-acoustic waves with Alfven waves in the solar wind, Lacombe and Mangeney came 
to the conclusion that "this non-linear process is not very efficient in the solar wind, so that 
Alfven waves can be considered as decoupled from compress ive waves in the major part of 



the m.h.d. spectral range" (see lLacombe fc Mangeney Ill980l . Abstract). Another argument 



against this scenario is that, if one of the modes is strongly damped, three- and/or four- wave 
resonant processes transform into the induced scattering of the involved waves by plasma 
particles (IBreizman. Zakharov fc Musher Ill973l ) and, hence, they do not lead to the reap- 
pearance of the damped wave mode. I n addition, ther e are observational arguments against 
the multi-wave scenario. It is known (IBiskamp 1120031 ) that in the slow solar wind the fluc- 
tuation amplitudes are smaller then in the fast solar wind. If nonlinear multi-wave resonant 
processes were responsible for the generation of compressive fluctuations, then one should 
expect a much higher level of density fluctuations in the fast wind than in the slow wind. 
However, observations s how that the slow solar wind is usually much mor e compressional 
then the fast solar wind ([Bruno fc Bavassano Ill99ll ; iBavassano et. al I I20Q4J ) . Consequently, 
the local generation of magneto-acoustic waves could hardly be explained by multi-wave 
resonant processes. 
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Hence, the gradual appearance of locally-produced magneto-acoustic waves throughout 
the solar wind seems to be an observationally confirmed fact, but still solicits for an adequate 
explanation. In the present paper, we argue that the velocity-shear-induced linear mode 
conversion of Alfven waves into fast and slow magneto-acoustic waves represents an efficient 
mechanism for the generation of the observed magneto-acoustic waves within the solar wind. 
It is also argued that the velocity shear could have an important contribution in the gradual 
"energization" of the solar wind via the generation of fast and slow magneto-acoustic MHD 
waves, that are eventually strongly damped by Landau damping. Our analysis shows that 
the linear mode conversion of Alfven waves into magneto-acoustic waves is especially efficient 
for the low-frequency Alfven waves from the en ergy containing interval. Afterwards, these 
perturbations can cascade to the smaller scales ([Montgomery et. al 1119871 ) by turbulence. 



It is well-known ( [Roberts et al. Ill992l ; iGoldstein fc Roberts Ill999l ) that velocity shear is 
one of the most important ingredients in the solar wind dynamics. Recent studies, based on 
in situ observations, showed that the fast/slow wind interface in the interplanetary space has 
two parts: a smoothly varying "boundary layer" flow that separates the fast wind from the 
coronal holes, and a shar per discontinuity between the slow and the intermediate solar wind 
( iSchwadron et al. II2005I ). A relatively high velocity shear was observed by Ulysses over its 
first orbit in the transition area between the fast and slow solar winds at 13° — 20° latitudes. 



The analysis of the data (jMcComas et al. Ill998l ) showed that the transition area between 
the two winds consisted of two regions, the first one with a width Al± rs 2 10 km and a 
velocity change AVi ~ 200 kms -1 , and the second one with Al 2 ~ 8 10 7 km and AV2 ~ 
100 kms -1 . However, during Ulysses' second orbit, the gl obal solar wind structure wa s 
remarkably different from that observed during its first orbit (jMcComas et al. 1 1200 ll . 120031 ). 
Overall, the solar wind was the slowest seen thus far in the satellite's ten-year journey ( 270 
km s -1 ). The wind was highly irregular, with less pronounced periodic stream interaction 
regions, more frequent coronal mass ejections, and only a single, short interval of fast solar 
wind. The complicated solar wind structure obviously was related with a higher complexity 
of the solar corona around the solar activity maximum, e.g. with the disappearance of large 
polar coronal holes and with the presence of smaller-scale coronal holes, frequent CMEs and 
coronal streamers. 

Originally, the idea of ( yelocity-)Shear-induced Wave Trans formations (SWTs) in the 
solar wind was proposed by iPoedts. Rogava fc Mahajan I (119981 ). However, the model of 
these authors was based on a phase-space analysis of the temporal evolution of individual 
fluctuation harmonics in the ideal (viscosity- and resistivity-less) limits, and it was only of 
a qualitative nature. Since then, considerable progress has been made in three important 
directions: (i) real (physical) space numerical simulations have been performed and it has 
been shown that SWTs occur in a well-pronounced way and they lead to easily recogniz- 
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able collective phenomena in MHD plasma flows (iBodo et al. Il200ll ); (ii) the role of the 
dissipation has been analyzed and the concept of shear-induced self-heating was introduced 
( IRogava 112004 ) and it was found that compressible wave modes (e.g., sound waves in hydro- 
dynamic flows and fast magneto-acoustic waves in MHD flows) grow non-exponentially and 
undergo subsequent viscous and/or resistive damping leading to the heating of the ambient 
flow by "inborn" waves; (iii) a noteworthy quantum-mechanical (QM) analogy has been dis- 
closed that helped to apply efficient mathematical tools from the scattering matrix theory 
of QM to the SWT studies which helped to give a quantitative rigor to this theory and to 
calculate directly the efficiency of different wave transformation channels by determining the 
corresponding transformation coefficients. 

The aim of the present paper is to apply the latter method to the study of MHD wave 
transformations in the solar wind and to give a fully quantified description of the coupling 
efficiency by calculating the corresponding transformation coefficients. The mathematical 
methods used in this paper are similar t o the ones that originally were d eveloped in th e 



1930s for quantum mechanical problems (jStuekelberg I Il932l ; IZener I Il932l ; lLandau 1 119321 ). 



More recently, similar asymptotic methods have been successfully applied to various other 



problems including the in t eraction of plasma waves in inh omogeneous media (ISwanson 111998 
Gogoberidze et al. 112004 ; IRogava k, Gogoberidze 1120051 ). 



The present paper has the following structure: the main mathematical consideration is 
presented in the next section. Different regimes of wave transformations, depending on the 
value of the plasma-/?, are studied in the third section. A brief discussion and the conclusions 
are given in the final section. 



2. Basic Formalism 

In this section, we give a brief s ynopsis of the basic shear flow model that was studied by 



Poedts. Rogava fc Mahajan I (119981 ). In this vein, we consider a plane-parallel, compressible, 
magnetized, unbounded shear flow with uniform values of the equilibrium plasma density 
(po) an d pressure (Pq). The equilibrium magnetic field B is considered to be uniform as 
well and is assumed to be directed parallel to the flow velocity which, in turn, is spatially 
inhomogeneous and oriented in the z-direction: 

U o = (0,0,Ar), (1) 

with the shear parameter A > being defined as a positive constant. 



The linearized ideal MHD equations governing the evolution of the perturbations of the 
plasma density (p'), pressure (p'), velocity field (u') and magnetic field (b') in this flow are 
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[with the notation T>t = d t + Uo • V]: 

V tP ' + p V • u' = 0, 

Po[V t u' + (u' ■ V)U ] = -W - -^B x (V x b') 
V t b' - (B • V)u' + B (V • u') = 0, 
V ■ b' = 0. 



(2) 

(3) 

(4) 
(5) 



The standard technique of the so-called "shearing-sheet approximation" (IGoldreich fc Lynden-Bell 



19651 ) implies that it is convenient to expand the perturbations as follows: 



*'(x, t) = $'(k, t) exp [i(k x (t)x + k y y + k z z)} 



(6) 



where the state vector <&' = [u', b', //, p'), while k x (t) = k x (0) — k z At, and k x (0), k y , and k z 
are the initial (t = 0) values of the components of the wave number vector k(t). Remark- 
ably, this ansatz reduces the mathematical part of the task to the solution of a closed set 
of ordinary differential equations (ODEs) in time: a solvable initial value problem. Note 
that k x (t) varies in time and this fact is usually referred to as the "drift" of the Spatial 



Fourier Harmonics (SFH) in the phase k-space (IChagelishvili. Rogava fc Tsiklauri 
Poedts. Rogava fc Mahajan Ill998l ). 



1996 



From the set of Eqs. (2)— (5), considering adiabatic perturbations p' = c s 2 p', where c s 
is the sound speed, and introducing the following non-dimensional parameters and variables 

R = A/{V A k z ), T EE k z V A t, 13 EE C 2 JVl K y EE k y /k Z) K x (t) EE k x /k z - Rt, K\t) = 

K 2 x (t) + K 2 y + 1, p(k, r) ee ip'(k, r)/p , b(k, r) ee ib'(k, r)/B , v(k, r) = u'(k, r)/V A , ^ = 
p + K x {r)b x + Kyby = p — b z , (Va = B n /\/4:7ipo is the A lfven speed), one can derive the 
following set of ODEs JPoedts. Rogava fc Mahaian Ill998h [F^ = d?F]: 



tfj {2) + C n il> + C 12 {r)b x + C 13 b y = 0, 

b x 2) + C 22 (r)b x + C 2 i(t)V + C 23 (r)b y = 0, 
bW +C 33 b y + C 31 ij + C 32 (r)b x = 0, 
where ||C|| is a 3 x 3 symmetric matrix defined as: 



P -PK x (r) 
-(3K x {r) l + (l + (3)K 2 x {r) 
-(3K y (l + (3)K x (T)K y 



-PK y 
[l + (3)K x {r)K y 
l + (l + P)K 2 



(7) 



(9) 



(10) 
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The total energy of the perturbations, in the non-dimensional form, can be defined as 
the sum of the kinetic, the magnetic and the compressional energies. It can be written down 
in the following way: 



E(k,r) 



|v| 2 + |b| 2 + /V 



From Eqs. (7)-(9) it is easy to see that in the absence of the shear in the flow (i.e. 
for R = 0, and thus A = 0), these equations describe indepen dent oscillations with the 
fundamental eigenfrequencies (IPoedts. Rogava fc Mahajan Ill998l ): 



1 ± 



4(3 



(1 + (3) 2 K 2 



Q 2 



(12a) 



(126) 



that can be easily identified as fast and slow magneto-acoustic waves (FMW and SMW) and 
Alfven waves (AW), respectivel y. For the eigenfunctions , ^j, corresponding to the above- 



defined fli eigenvalues, we find (IGogoberidze et al. 112004 ) 



_ (n 2 s -/3)p + n 2 s (K x b x + K y b y ) 

(5K\p + (Q 2 - (3K 2 )(K x b x + Kyby) 
K ± ^(n 2 -(3) 2 + (3 2 K 2 ± 



Ky h x K x by 



(13a) 

(136) 
(13c) 



where K 2 = K 2 + K 2 



Notice than, when the flows are only weakly sheared (i.e. when the coefficients 

in Eqs. (7)-(9) vary only slowly or adiabatically. This implies that the expressions for the 
fundamental eigenfrequencies, given by Eq. (12), and the corresponding eigenf unctions, given 
by Eq. (13), are still useful for a qualitative description o f the shear-induced dynamics of the 
wave modes and their coupling /conver sion properties (jChagelishvili. Rogava fc Tsiklauri 
19961 : IPoedts. Rogava fc Mahaian Ill998h . 



3. Transformation regimes: the transition matrix method 



As will be shown in the next section, the characteristic values of the normalized velocity 
shear rate in the solar wind plasma always satisfy the condition R <^ 1. This necessarily 
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implies that the coefficients in Eqs. (7)-(9) are only slowly varying functions of r and, 
therefore, the adiabatic (WKB) approximation holds everywhere except in the immediate 
vicinity of the turning [Qi(r t ) = 0] and the resonant [fii(r r ) = f2 2 (T r )] points. Using Eq. (12) 
one can evaluate that the condition: 

tu < n?, (14) 

is satisfied for all the MHD wave modes at any moment of time, or equivalently, none of the 
turning points are located near the real r-a xis. From a physical point of view, this means that 



there are no (over-)reflection phenomena (IGogoberidze et al. 



2004 ) and resonant coupling 



can only occur between different waves modes with the same sign of the phase velocity. 

A careful analysis of the system yiel ds that the resonant coupling tak es place in the 
vicinity of the point r* where K x (t*) = (IPoedts. Rogava fc Mahajan Ill998l ). According to 
the general theory o f such systems, the timesc ale of the resonant coupling At is of the order 
of At ~ i2~ n /( n + 1 ) ( Gogoberidze et al. 112004 ). where n indicates the order of the resonant 
point 0, and, therefore, the evolution of the waves is adiabatic when 



\K x (t)\ >i? 1/(n+1) . 



(15) 



If this condition is satisfied, the temporal evolution of the waves is described by the 
standard WKB solutions: 



3 ±i/f2i(r)dT 



(16) 



where the D i denote the WKB amplitudes of the wave modes with positive and negative 
phase velocity along the z-axis, respectively. All the physical quantities can be readily found 
by combining the Eqs. (7)-(9). The energies of the involved wave modes satisfy the standard 
adiabatic evolution condition (IPoedts. Rogava fc Mahajan Ill998l ): 



^ = ^(r)(|^ 



(17) 



From this equation it follows that l-D^] 2 can be interpreted as the number of 'wave particles' 
(the so-called 'plasmons'), in analogy with quantum mechanics. 

Let us assume that initially K x (0) R 1 ^ n+1 \ Due to the linear drift in the k-space, 
K t (t ) decreases and when the mode enters the "degeneracy area" (IPoedts. Rogava &: Mahajan 
19981 ). the mode dynamics becomes non-adiabatic due to the resonant coupling between the 
modes. Afterwards, when K x {r) <C — i? 1//( ^ n+1 \ the evolution becomes adiabatic again. De- 
noting the WKB amplitudes of the wave modes before and after the coupling region by 



1 The resonant point is said to be of the order n when (Of 



(r — T r )™/ 2 in the neighborhood of r r 
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D iB and D i A , res pectively, and by em ploying the formal analogy with the ^-matrix of the 

Z A with Df B via the so-called 



scattering theory ( Kopaleishvili I Il995l ). one can connect 
transition matrix: 



D 



TD 



B- 



D 



T*D 



B ' 



where T and its Hermitian conjugated matrix T* are 3x3 matrices. Note, that none of the 
turning points are located near the real r-axis and, therefore, there is n o transition between 



the m odes with opposite signs of the phase velocity along the z-axis ( jGogoberidze et al. 
2004h . 



Notice that all the coefficients in the governing equations are real. Moreover, the m atrix 
given by Eq. (10) is symmetric. As a consequence (jKopaleishvili Ill995l ; iFedoriuk Ill983l ). the 
transition matrix T is unitary 0, and 



(19) 



Generally speaking, this equation represents the conservation of the wave action. When 
R <C 1, it tra nscribes into the energy c onservation throughout the resonant coupling of the 
wave modes (jGogoberidze et al. II2004I ). Energetically this means that a transformed wave 
mode is generated solely on the expense of the energy of the incident wave mode. 

The crucial physical importance of this matrix follows from the fact that the value of 
the quantity \Tij\ 2 represents a part of the energy transformed via the resonant coupling of 
the modes. That is why the absolute values of the transition matrix components |T^| are 
called the 'transformation coefficients' of the corresponding wave modes. For the resonant 
interaction of two wave modes, e.g. i and j, the unitarity of the matrix T guarantees an 
important property of the transition matrix, viz. Tij = Tji. This symmetry property holds 
for the resonant interaction of two wave modes only. If, in the same time interval, there is 
an effective coupling of more then two wave modes, then the symmetry property fails. 

I n earlier studies, a lot of attention was usu ally paid to the resonant points of the first 
order (ILandau &: Lifschitz 1119771 ; IFedoriuk Ill983l ). In this case, only the dispersion equations 
of the waves are needed to derive the transformation coefficients with accuracy 0{R 1 ^ 2 ). For 
the second and/or higher order resonant points, analytical expressions for transformation 
coefficients can be derived only in the case of weak interactions (T^ < 1, i / j). As a 
matter of fact, if the resonant points are not close to the real r axis, in the sense that 



'13 



Im / (fij — Vtj)dr 



TO 



> 1, 



(20) 



2 A matrix U is unitary if its conjugate transpose, U H is equal to the inverse matrix U 1 . 
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then the transformation coefficient is just equal to: 

Tij » - exp (-0^) . (21) 

Here, and hereafter, the signs of the absolute magnitude for transformation coefficients are 
omitted, i.e., from now on the notation means |Ty|. 

In the context of the linear problem that we are studying in this paper, viz. the velocity- 
shear-induced coupling of MHD waves, all the resonant points are of the second order, as 
will be shown later. 

In the following subsections, we will study in detail the coupling of MHD wave modes 
(AW, FMW, and SMW) for different regimes of the plasma-/^ that are interesting in the solar 
wind context. For the earlier , detai led, but only qualitative analysis of the same regimes see 



Rogava. Poedts & Mahaian I (l2000h . 



3.1. The case (3 < 1 

It is well-known that the magnetic field dominates the plasma (/3 <C 1) throughout most 
of the solar corona, especially at lower altitudes and within the so-called 'active' regions, as 
well as in the innermost regions of the solar wind. For the MHD waves in this regime, 
the frequency of the SMW is far smaller then the frequencies of the FMW and the AW. 
Therefore, the coupling of the SMW with the other two MHD modes is exponentially small 
with respect to the large parameter 1/R and can be neglected. Consequently, in the set 
of the governing equations Eqs. (7)-(9), the equation for the variable if) decouples from the 
other equations, and the equations for b x and b y describe the evolution of the coupled AW 
and FMW: 

b x + [1 + K 2 x (t)} b x = -K x {r)K y b y , (22) 
b y + [1 + K 2 y ] b y = -K x (r)K y b x . (23) 

The normalized frequencies of the coupled wave modes are: f^(r) = l+Ky+K^(r), fl\ = 
1. Therefore, there are two, complex conjugated, second order resonant points: 

K x {r fA ) = iK y , K x (r* fA ) = -iK y . (24) 

The condition (15) then implies that the evolution of the waves is adiabatic if 

6=\K%\/R^>1, (25) 
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and if this condition is satisfied, Eq. (21) yields: 



T 



7T 



fA 



2 ex n _ y 



(26) 



An analytical expression for t he transformation coefficie nts can also be derived in the 
opposite limit, 5< 1. In this case (IGogoberidze et al. 1120041 ): 



T 



2 2/3 



7T 



fA 



3vs r (i) 



6 1 



r(|) 



2 7/4 3 l/3 r (|) 



S' 1 



(27) 



Results of the numerical solution of the initial set of equations (7)-(9) (solid line) as 
well as of the analytical expressions given by Eqs. (26) (dash-dotted line) and (27) (dashed 
line) are presented in Fig. 1. This plot shows that the transformation coefficient attains its 
maximal value (Tj A ) max = 1/2 at 5 cr « 0.89. 

In the case of first-order resonant points, the transformatio n coefficient monotonically 
tends to unity when the resonant points tend to the real axis (ILandau &: Lifschitz 1 119771 ; 
Fedoriuk I Il983l ) and, therefore, the total conversion of one wave mode into another wave 
mode is possible. Here, in the case of second-order resonant points, existing for the FMW- 
AW coupling in a low-/? plasma, (TfA 2 ) ma x = 1/2. This leads to the following important 
astrophysical conclusion: in those regions of the solar wind where (5 -C 1, even under the 
most favorable conditions, at most half of the energy of the AWs can be transformed into 
FMWs and vice versa! 

Yet another important aspect of the coupling in this case is that, if the resonant point 
tends to the real r axis (i.e., K y — > 0), the transformation coefficient tends to zero. 



3.2. The case (3 > 1 

In the high-/? case, Q s , Qa *C fi/. Therefore, only the coupling between the AW and the 
SMW can be important in this regime. From the expressions for the fundamental frequencies 
[see Eq. (12)] one can easily deduce that, for the AW-SMW coupling, there are two complex 
conjugated, second-order resonant points, that are given by Eq. (24). The condition (20) in 
this case has the following form: 



1 

~f3R 



arcsinh (K y 



> l. 



(2? 



12 



If this condition is satisfied, and in addition K y <C 1, we can derive from Eq. (21) the 
following formula for the transformation coefficient: 



Unlike in the low-/? case, now there is no unique parameter, like the parameter 5, which 
would give a complete description of the transformation process. Instead, in this limit, there 
are two "governing" parameters, viz. K y and j3R. 

Let us first consider the case f3R <C 1. In this case, Eq. (28) reduces to 8\ = \K y {/(fiR) 1 / 3 ^> 
1, and the properties of the wave transformation are essentially the same as in the case of 
the transformation of the AW to the FMW and vice versa. As a matter of fact, if Si <C 1, 
then the leading term of the asymptotic expressions of the transformation coefficient is given 
by Eq. (27), with 5 replaced by Si. At f3R <C 1, the transformation coefficient reaches its 
maximum (T^ A ) max = 1/2 at 5f that coincides with 5 cr . 

The dependence of the transformation coefficient on K y , on the other hand, provides 
some very interesting and new details, compared with the low-/? case. On Fig. 2, we display 
the numerical solution of the initial set of equations (7)-(9) for the cases f3R = 0.025 and 
f3R = 1. This numerical inspection shows a remarkable fact: when j3R is not small the 
properties of the transformation process are quite different (see Fig. 2). As a matter of fact, 
when K y 1 , it turns out that the transformation coefficient T s a does not depend on f3R at 
all and is given by the formula T s a ~ 2.05K y . Besides, as it can be seen clearly from Fig. 2, 
{T S A)max = 1, i-e. now, unlike the low-/? situation, a total transformation of one wave mode 
into another wave mode is possible! 



The (3 ~ 1 case is the most interesting but also the most complicated case of velocity- 
shear-induced MHD wave transformations. The frequencies of all three modes are of the 
same order in this case and a simplification or reduction of the set of Eqs. (7)-(9) is not 
possible. However, the analysis of Eq. (12) enables us to see that there is a pair of complex 
conjugated first-order resonant points: 




(29) 



3.3. The (3 ~ 1 case 




(30) 



and another pair of complex conjugated resonant points of the second order: 



K x (t 3A ) = ±iK v . 



(31) 
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The essential complexity of this case stems from the fact that the transition matrix 
method does not allow one to derive analytic expressions for the transformation coefficients 
when more than two wave modes are effectively coup l ed. H owever, the numerical study of 
the problem shows (see for details iGogoberidze et al. I (120041 )) that the qualitative character 
of the wave transformation processes in this case is mainly similar to the cases described in 
the previous subsections. 



4. Discussion and conclusions 



Several years ago. lPoedts. Rogava &: Mahajan I (119981 ) argued that velocity-shear- induced 



MHD wave transformations could be an important ingredient of the wave d ynamics in the so 



lar w i nd, contributing to its acceleration a nd heating. It was further argued (jPoedts. Rogava fc Mahajan 



19991 ; iRogava. Poedts fc Mahajan 1120001 ) that SWTs can contribute to the transmission of 



the waves from the chromosphere, through the transition layer to the corona and further 



(Rogava 


2004; 


Shereelashvili, Roeava. & Poedts 


2005: 


Shereelashvili. Poedts & Patarava 


2006 


), via the agency of mutually coupled wave modes, might be one of the reasons for the 



coronal heating (both "basal" and "extended" ) and the acceleration of the solar wind. 

Up to now, the scenario of the linear coupling and the mutual transformations of the 
MHD wave modes in the solar wind was outlined only qualitatively. Clearly, this merely 
qualitative picture lacked quantitative strength and rigor, since it was not clear how efficient 
the SWTs could be in the solar wind context. In this paper, we solved this problem and 
we provided a systematic, quantitative description of SWTs in terms of the recently devel- 
oped transition matrix method (originally developed in the 1930s for quantum mechanical 
applications) and transformation coefficients. 

In particular, we have demonstrated that the dynamics of the MHD wave conversions 
is determined by three key parameters: the flow velocity shearing rate R, the plasma f3 and 
the ratio K y = k y /k z . Since the aim of this paper is to verify whether the SWTs are well- 
pronounced in the solar wind, it is necessary to analyze whether these parameters within the 
solar wind really have values which favor efficient mode transformation processes or not. 



For the dimensionless shearing parameter R we have 

R AV1 



Al tt a 



(32) 



Using for AV and Al the values presented in the Introduction, we found that even for very 
low frequency Alfven waves (with Q a > 10 -4 s _1 ), the dimensionless shear parameter satisfies 
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the condition R -C 1, for both regions of the transition area between the fast and slow solar 
winds. Bearing also in mind that the plasma (3 ~ 1 in the solar wind, we obtain that Rf3 <C 1, 
and according to the considerations in the previous section, the coupling of the AWs with 
both the fast and slow magneto-acoustic waves has the same qualitative character. As a 
matter of fact, Eq. (25) gives us the necessary condition for an effective coupling 

K„ < R l l\ 



(33) 



showing that an effective coupling between the Alfven mode and the fast and slow magneto- 
acoustic modes takes place for the low frequency Alfven waves from the energy containing 
interval. Indeed, using for the frequency of the AWs, Q a = 2tt /T a k6x 10~ 4 s _1 , and using 
for AV and Al the values presented in the Introduction, we found that Ri « 1.6 x 10~ 2 
and i?2~2x 10~ 3 for the two regions of transition between the fast and slow solar winds. 
Equation (33) then yields for the perpendicular wave numbers of the waves, k y < 0.25k z 
for the first transition region and k y < 0.l3k z for second transition region, respectively. 
Noting that in the e nergy containing range, the turbulent fluctuations are nearly isotropic 
( iBigazzi et. alll2006l ). we can expect that for the typical turbulent fluctuations in this range 



parallel and perpendicular wave numbers of the same order [k y ~ k z ). Consequently, the 
linear transformations of Alfven waves to fast and slow magneto-acoustic waves are quite 
likely to take place in this frequency band. 

On the other hand, the linear transformation of Alfven waves to magneto-acoustic waves 
seems to be inefficient in the high frequency band, including the inertial interval of solar 
wind fluctuations. As a matter of fact, Eqs. (32) and (33) yield that the modes which 
could be effectively transformed should satisfy k y <C k z . However, theoretical research 
of MHD turbulence as well as numerical simulations and analysis of the solar wind data 



(Shebal in et al. 



1983 



Golderich fc Sridhar 



1995 



Cho fc Vishniac 2000: Miiller et al. 2003 



Biskamp II2003I : lOughton fc Matthaeus II2OO5I " iGogoberidze 1120061 ) all indicate that the energy 



cascade in the inertial interval proceeds much more effectively in the direction perpendicular 
to the mean magnetic field and, therefore, for high-frequency modes one usually has k y 3> k z . 
Consequently, the condition (33) could be hardly fulfilled for typical parameters of high- 
frequency Alfven waves in the solar wind. 

Although the performed analysis allows us to make quite definite predictions about the 
efficiency of linear conversion of AWs into FMWs and SMWs, direct quantitative comparison 
with the observations is difficult since the compressive par t of the velocity field can not be 
isolated by single spacecraft observations (IBiskamp 1120031 ). Consequently, it is difficult to 
associate the observed density fluctuations to specific mode characteristics, such as fast or 
slow magneto-acoustic modes. 



Nevertheless, on the qualitative level, the predictions of the study presented here are in 
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good agreement with the observations. As a matter of fact, (i) in the regions with strong 
velocity shear , Alfvenic correlations are r educed and a high level of density perturbations 
are observed ([Goldstein fc Roberts Ill999l ); (ii) the spectrum of density fluctuations in the 
i nertial interval of t he solar wind fluctuations follows the Kol mogorov scaling p' ~ fc -5 / 3 



(IMarsch fc Tu II19901 ). This agrees with the model developed by Montgomery et. al I (119871 ). 
where it has been shown that in the weakly compressible limit, the density perturbations 
behave like a passive scalar, i.e., they follow the scaling of the magnetic field perturbations. 
The findings of the present paper supplement this model and offer an efficient mechanism 
for the generation of the density perturbations in the energy containing interval, which 
afterwards are cascaded to the smaller scales. 

Another interesting and not completely understood problem is the spatial aspect of the 
velocity-shear-induced wave transformations in the solar wind. If the largest values of the 
flow shearing rates appear in the relatively narrow "transition zones" between the slow and 
the fast solar wind, one may wonder whether the velocity-shear-induced mode conversions are 
confined to these regions or whether they tend to occupy a wider volume of the wind plasma!? 
Evidently only real-space numerical simulations, taking into account physical characteristics 
of the solar wind, may answer this ques tion. Direct numer ical simulations, performed in the 
general context of SWTs, have shown (IBodo et al. Il200ll ) that when the Alfven waves get 
transformed into the fast magneto-sonic waves, the latter tend to leave their "birthplace" 
and propagate quite rapidly due to their high velocities. We suppose that the observational 
signature of the waves generated by means of this mechanism would be their appearance 
throughout both the sheared and the non-sheared wind volume with the maximal probability 
of their detection in the vicinity of the strongly sheared transition regions. 

In t he near future, we would l ike to take into account the effect of the kinematic com- 
plexity (IMahajan fc Rogava I fl99~9l ) and verify whether the transition matrix method can 
lead to a quantitative description of the SWTs not only for the solar wind, but also for solar 
jet-like flows with a more complicated geometry and complicated kinematics. 
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Figure captions 

FIG. 1: Transformation coefficient Tf A vs 5. Dash-dotted line and dashed line represent 
analytical expressions (26) and (27), respectively. The solid line is obtained by the numerical 
solution of Eqs. (7)-(9). 

FIG. 2: Transformation coefficient T s a vs K y for R = 0.005 and /3 — 5 and (3 = 200. 



-20- 




Fig. 2.- 



